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Abstract: We reformulate perturbation theory for neutrino oscillations in matter with 
an expansion parameter related to the ratio of the solar to the atmospheric Am^ scales. 
Unlike previous works, we use a renormalized basis in which certain hrst-order effects are 
taken into account in the zeroth-order Hamiltonian. We show that the new framework has 
an exceptional feature that leads to the neutrino oscillation probability in matter with the 
same structure as in vacuum to hrst order in the expansion parameter. It facilitates imme¬ 
diate physical interpretation of the formulas, and makes the expressions for the neutrino 
oscillation probabilities extremely simple and compact. We hnd, for example, that the i^e 
disappearance probability at this order is of a simple two-flavor form with an appropriately 
identified mixing angle and Am^. More generally, all the oscillation probabilities can be 
written in the universal form with the channel-discrimination coefficient of 0, ±1 or simple 
functions of 023 • Despite their simple forms they include all order effects of 0i3 and all order 
effects of the matter potential, to hrst order in our expansion parameter. 
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1 Introduction 

Neutrino oscillation based on the standard three-flavor scheme provides the best possi¬ 


ble theoretical framework available to date to describe most of the experimental results 


obtained in the atmospheric, solar, reactor, and the accelerator neutrino experiments. 
Although numerically calculated neutrino oscillation probabilities suffice to analyze exper¬ 
imental data, understanding of the framework, in particular the one in matter [1], has not 
yet reached a sufficient level, in our opinion. Under the assumption of uniform matter den¬ 
sity distribution, the exact expressions of the eigenvalues, mixing angles, and the oscillation 
probabilities in matter have been obtained [2-4]. Yet, the results for these quantities are 
generally too complicated to facilitate understanding of the structure of the three flavor 
neutrino oscillations in matter primarily due to the complexities of solving the cubic eigen¬ 
value characteristic equation. For a recent comprehensive treatment of neutrino oscillation 
in the matter, see ref. [5]. 

Analytic approaches to the neutrino oscillation phenomenon, so far, are mostly based 
on variety of perturbative frameworks. If the matter effect is small one can treat it as a 
small perturbation [6]. In the environments in which the matter effect is comparable to 
the vacuum mixing effect, the only available small expansion parameter known to us is the 
ratio of the solar-scale Atuq to the atmospheric-scale Am^, Atuq/Atu^ ~ 0.03. sin^ia 
has been often used as an expansion parameter (there are enormous number of references, 
see e.g., [7]), but it is now known that its value is not so small, sin0i3 ~ 0.15, which is of 
the order of y^Am^^/Am^. Moreover, expansion around sin0i3 = 0 misses the physics of 
the resonance which exists at an energy around ~ 10 GeV for earth densities. Therefore, 
it appears that the suitable perturbative framework is the one with the unique expansion 
parameter AmQ/Am^. This framework was indeed examined in the past, to our knowledge 
in refs. [7-10]. 

In this paper, we present a new framework of perturbative treatment of neutrino 
oscillation in matter. We follow the reasoning stated above which led to identification of 
the unique expansion parameter e ~ Am-Q/Am^. But, unlike the preceding works, we use a 
“renormalized basis” as the basis of perturbation. That is, we absorb certain terms of order 
e to our “zeroth-order” Hamiltonian around which we perturb. Or, in other word, we take 
the zeroth-order eigenvalues in matter such that it matches the exact eigenvalues to order 
e, see section 5. We will show that use of the renormalized basis makes the structure of the 
perturbation theory exceptionally transparent, as we will explain in the next section. It 
allows us to obtain simple, elegant and compact expressions for the oscillation probabilities, 
which have a universal form even to first order in our expansion parameter. 

For example, Ve survival probability takes the form to order e as 



( 1 . 1 ) 
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where (j) is 6*13 in matter, and X± denote the eigenvalues of the states which participate 
the 1-3 level crossing. Despite its extremely simple form, P{i'e ^e) in (1-1) takes into 
account all order effects of @13 and the matter potential. Since we will only consider terms 
up to order e, in this paper, our results here are not applicable to the region near the solar 
MSW resonance [1, 11]. Our perturbative framework will be called as the “renormalized 
helio-perturbation theory” in the rest of this paper. 

The section plan of this paper is somewhat unusual: in the next section 2 we summarize 
the key features of the perturbative framework we develop in this paper. Then, in section 3 
we describe the principle results of this work including the oscillation probabilities for all 
channels in matter. This section does not describe the derivations but provides a self 
contained summary of the results of this paper. Following this section, see section 4, we 
present a systematic exposition of our perturbative framework and how the results of the 
section 3 are derived. In the appendices A, B, and C we present, respectively, calculational 
details of the S matrix, the results of oscillation probabilities to order e in the standardized 
form, and useful relationships to verify the equivalence of various expressions. 


2 Structure of the neutrino oscillation probabilities in vacuum and in 
uniform matter 


We start by giving a precise definition of what we mean by “simplicity and compactness” 
of the expressions for the neutrino oscillation probabilities. Then, we explain the need for 
reformulating perturbative treatment, followed by outlining “what is new” in this paper. 


2.1 Simplicity and Compactness 

In vacuum and in matter with constant density, it is well known that the neutrino oscillation 
probabilities for vp —?• Va for three-flavor mixing {i,j = 1,2,3) can be written in the 
following form^ 




Vn) = 




i=l 


sin 

j>i 


2 i^j Ai)L 


AE 


I CT FT/ T/* TA* TA 1 • (^3 - A2)T . (A2-Ai)L . (Ai - A3)L 
+ o lm[VaiVpiV^2^l32\ sin- jrpp -sm- — -sm- — -. ( 2 . 1 ) 


AE 


AE 


AE 


The flavor mixing matrix elements, VAi, relate the flavor states, i/q, to the eigenstates of 
the Hamiltonian, t'j, with eigenvalues Xij^lE by = VaiVi- Equation (2.1) applies both in 
vacuum and in uniform matter since the matrix V diagonalizes the full Hamiltonian, which 
includes the Wolfenstein matter potential, a [Ij. In vacuum. A* = m? where denotes 
the mass of the z-th neutrino state and V = Umns- 

Notice that the identical form of the oscillation probabilities in vacuum and in matter, 
eq. ( 2 . 1 ), apart from replacement of Amfj by (Aj — Xj), implies that physical interpretation 
of the formulas in matter is as transparent as in vacuum. 

^ We have used for the last term in ( 2 . 1 ) W2V'/32] = Iui[Vij2F^*2W3E33] = W3V83] 

which follows from unitarity, and the identity (C. 2 ). 
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What is less well appreciated is that the expressions of the oscillation probabilities in 
(2.1) are maximally simple and compact. That is, they contain 5 (including the constant 
term) functions oiLjE which are linearly independent. This property can be easily verified 
by showing the Wronskian is nonvanishing, which implies that none of the functions can be 
written in terms of the other functions for all L/E. Hence, in a true mathematical sense, 
eq. (2.1) and it’s equivalents give the simplest and most compact form for the complete set 
of three flavor oscillation probabilities in vacuum and in uniform matter.^ 

A special feature of the oscillation probabilities which is also worth noting is that 
each term in eq. (2.1) factorizes into the characteristic sin[(Aj — \i)L/AE] factor and the 
products of the V matrix elements which control the amplitude of the oscillation. Both the 
eigenvalues. A* and the matrix elements of V are independent of the baseline, L, but are 
functions of the mixing angles, 0’s, the Am|j, and the product of the energy of neutrino 
times the matter density via the matter potential.^ CP and T violation is described by the 
last term in (2.1), which has the universal, channel-independent form in the three neutrino 
mixing. 

2.2 Raison d’ etre and the requirements for perturbative treatment 

With the simplest form of the oscillation probability eq. (2.1) and by knowing both the 
exact form of the eigenvalues, Aj/2H, and the elements of the V matrix, see [3], one might 
expect we have all that is needed for theoretical discussions. Unfortunately, the analytic 
expressions for these eigenvalues, Aj/2U, as well as the V matrix elements are notoriously 
complex and give no analytic insight into the oscillation physics in uniform matter. This 
is true even when one of sin0i3, sin0i2 or a is set equal to zero. In any one of 

these limits, the characteristic equation for the eigenvalues factorizes. But, the form of the 
general solution does not simplify trivially to yield the correct eigenvalues, even though 
it must. The structure of the general solutions of the cubic characteristic equation is the 
root cause of this rather unusual behaviour. Hence, there is a need for a reformulated 
perturbative framework so that we can obtain approximate but much simpler expressions 
for the eigenvalues, Aj, and the mixing matrix elements, Vai, which provide the necessary 
physics insight. 

In this paper, we formulate perturbation theory by which we can calculate the eigen¬ 
values Aj and the elements of V matrix as a simple power series of the small expansion 
parameter, a renormalized Attiq/A m'^. At the same time the structure-revealing form 
of the oscillation probabilities (2.1) is kept intact. While the existing frameworks do not 
satisfy the latter requirement, the key to the success in our case is due to the correct 
decomposition of the Hamiltonian into the unperturbed and the perturbed terms. Use of 

^ There are equivalent ways to write this set of oscillation probabilities with 5 independent L/E functions, 
e.g one could use the trigonometric identity 2 sin^ x = 1 — cos 2x to replace the sin^ x. Any other way of 
writing these oscillation probabilities will have 5 or more linearly independent functions, e.g., using the 
identity (C.2) increases the number of independent functions by 2. 

® a oc pE with p being the matter density. The V matrix elements are determined solely by the elements 
of the Hamiltonian matrix. It is easy to observe this feature by using the method developed in ref. [4], 
which is valid for nonuniform matter density as well. 
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the renormalized zeroth-order basis, allows us to correctly determine the eigenvalues to the 
appropriate order in our expansion parameter. The resultant renormalized eigenvalues and 
the mixing parameters effectively absorb the additional terms that arise in the conventional 
perturbative frameworks, as extra functions to the minimally required 5 linearly indepen¬ 
dent functions. This occurs automatically and we believe that such a framework has never 
before been formulated.^ Because of the structural simplicity of the LjE dependence, the 
expressions of the oscillation probabilities calculated by our method are extremely simple 
and compact, as will be fully demonstrated in the next section. 

3 Results of our perturbative expansion for all oscillation probabilities 

In this section, we describe the main results of this paper without derivations and with 
only minimal discussion. In later sections we provide the derivation and more in depth dis¬ 
cussions. We start with the approximate eigenvalues of the Hamiltonian, the approximate 
neutrino mixing matrix and then give the oscillation probabilities for all channels to first 
order in the expansion parameter, e, see eq. (3.3) for the precise definition. 

3.1 Mass squared eigenvalues in matter 

In vacuum the three eigenvalues of the full Hamiltonian which governs the neutrino oscilla¬ 
tion is given in the form m? /2E, where rrii is the mass of i-th mass eigenstate of neutrinos, 
i = 1, 2, 3. Similarly, in matter we write the three eigenvalues as 


2E 


where the state label runs over i = — ,0,-|- for the approximate Hamiltonian of three 
flavor mixing system. To treat the normal and the inverted mass orderings (NO and 10 


respectively) in a unified way, we define the eigenvalues as follows^ 



Ao = c?2 e Am^en, 


(3.1) 



'ren'’12' 


^ See, however, the clarifying remark at the end of subsection 3.3.6. 

® We note that the eigenvalues in (3.1) above appear in ref. [5]. See section 4.1 for the derivation and a 


comment on the treatment in [5]. 
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In eq. (3.1), the renomalized Am^ = Am^g^, the expansion parameter e, and the Wolfen- 
stein matter potential [ 1 ], a, are defined as follows:® 

Am^en = Amii - s^sAm^i, 

Amni 


and a = iV^GpNeE ^ 1.52 x 10“^ ( A (^ ) eV^ 

\g.cm \GeV / 

This choice of Am^^^ is crucial to the compact formulas for the oscillation probabilities that 
will be given in this paper. Note also that the sign of Am^g^ signals the mass ordering, both 
Am^gj^ and e are positive (negative) for NO (10). However, for both orderings eAm^g^ = 
Am^i > 0 , as required by nature. Notice that Aq is the same for the both mass orderings, 
and when we switch from NO to 10 we also switch the sign in front of the square root 
in eq. (3.1). The nicest feature of the sign choice is that the oscillation probability has 
a unihed expression and the solar resonance is in level crossing for the both mass 

orderings. 

Am^en is equal to the effective atmospheric Am? measured in a electron (anti-) neutrino 
disappearance experiment in vacuum, Am^^ = 0 ^ 2 Aml^ -|- sf 2^^32 [l^]- This quantity 
is identical to Am^^ recently measured by the reactor @13 experiment [13] up to effects 
of 0{Am2i/Am'^i)'^. Whether the coincidence between Am^^ and reflects a deep 

aspect of neutrino oscillation or not will be judged depending upon what happens at second 
order in e. This point as well as the relevance of the other effective Am?^^ [12], equivalent 
of ArUgg, will be discussed in depth in a forthcoming communication. 


(3.2) 

(3.3) 

(3.4) 


3.2 The mixing angle @13 and mixing matrix in matter 

We use the angle cj) to represent the mixing angle @13 in matter. With the definitions of 
the eigenvalues (3.1), the following mass-ordering independent expressions for cosine and 
sine 2(f) (see section 4.2) are given by 


cos 2(j) 


sin 2(j) 


Am^gjj cos 2013 — CL 
X+-X- ■ 

Am^gn sin 20 i 3 
A+ - A_ ■ 


(3.5) 


It is easy to show that (p goes from 0 —?■ tt/ 2 as a goes from — 00 to -|- 00 for the NO 
and as a goes from -|- 00 to — 00 for the 10. In vacuum (a = 0), cp = 0i3 and (p = tt/A at 
the atmospheric resonance, when a = Am^g^ cos 20 i 3 , for both mass orderings. 

® The following notation is nsed throughout: Amfj = mf — m?, Sij = svnOij and dj = cos^ij where 
Oij are the standard neutrino mixing angles and Gf is the Fermi constant, W is the number density of 
electrons, E is the energy of the neutrino, the electron fraction and p is the density of matter. 
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The mixing matrix in matter, V, relates the flavor eigenstates, i^e, to the 

perturbatively defined matter mass eigenstates, as follows (see section 4.4): 


\^r ) 



(3.6) 


where the matrix V is unitary. It is convenient to split V into a zeroth order term, 
and a first order term, in our e expansion. 


= +eyW, 


(3.7) 


where the zeroth order matrix is given by 




—'S0C23 


0 S(p 

C23 C0S23e*'^ 

- 5236 “**^ C^C2‘i 


(3.8) 


with 5 being the CP violating phase, whereas the first order correction is given by 


y(i) 


C 12 S 12 < 


f a 

Va- -Aoy 


{ '^(0-613) A 

VA+-Aoy 


0 -C0 0 

C23 s</.S23e*^ 0 

-5236“*'^ S0C23 0 


'o 

Sip 

0 

' 

0 

-£05236*"^ 

C23 


0 

— C0C23 




(3.9) 


A factorized form of the V matrix is given in eq. (4.25). 

As an outcome of the consistent perturbative treatment the total V matrix given by 
(3.7) with (3.8) and (3.9) must be unitary to order e. In fact it is, since the following two 
conditions are satisfied 


y(o)(y(o))t = 1 and p(°)(p(^))'f + p(^)(I/(°))'f = 0. (3.10) 


Of course, none of what follows is self consistent without unitarity here. 

With the matter eigenvalues, A’s , definite by eq. (3.1) and the matter mixing matrix, V, 
given by eq. (3.7), simple and compact expressions can be easily derived for the oscillation 
probabilities in matter for all channels, to leading order in e, as will be shown in the next 
section. 


3.3 Compact formulas for the oscillation probabilities in matter 

In this section we start by presenting the shortest path to the oscillation probabilities of 
the Ve-TeXsXed channel by using the eigenvalues, and mixing matrix, V, given in the 
previous section to order e. Then, we derive the universal expression of the oscillation 
probabilities which is applicable to all channels. 
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3.3.1 Ve ^ disappearance channel 

The derivation of the Ve survival oscillation probability, P{i'e —^ i^e), in our renormalized 
helio-perturbation theory is extremely simple. Starting from the general expression 

P(/^e ^ I'e) ^ 1 - 4H4+|^|V;_|^sm^ 

where, L, is the baseline. Now |14oP = so we obtain to order e eq. (1.1), or 


P{ve Ue) = I- sin^ 26 »i 3 ^ sin^ (3.11) 

where |A+ — A_| = ^ (Am^g^ — a)^ + Is^gaAm^gn from eq. (3.1). 

Notice that the formula in eq. (3.11) takes into account the matter effect as well as the 
effect of si 3 to all orders. Nonetheless, it keeps an exceptional simplicity, an effective two- 
flavor form in matter which consists of single term with the unique eigenvalue difference 
A+ — A_, the feature we believe to be unique in the market. The feature stems from the 
fact that there is no Uf, component at zeroth order in e in the “0” state in matter. It is 
expressed in the zero in the 14o element of the zeroth-order V matrix as in (3.8), see also 
section 4.4. 


3.3.2 i^e —)• and Vg —)> Vr appearance channels 

Now, we discuss the appearance channels t'e and Ue ^r- We describe here the 

simplest way to derive the formulas for the oscillation probabilities starting from the V 
matrix by using unitarity. The oscillation probability P{i'e —^ i^fi) can be computed as 


Pil^e 


V„) = 




(3.12) 


We use unitarity relation V^-V*_ -|- l^oKo + '^^^+^e+ = 0 to eliminate the V^-V*_ term in 
(3.12). Then, we obtain 

P{ye ^ = 4|y^+lA;+sinA+_e-'^+“ - sin A_op 

= 4|y^+|2|y,+ |2sin2A+_ 

-STZ{V^+V*j^V*QVeo) sin A+_ sin A_o cos A+o 
-^I{V^+V:+V;^Veo) sin A+_ sin A_o sin A+o 
+MVM^\Veo\^ Sin^ A_o (3.13) 
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where the common shorthand notation for the kinematic phase Ajj = (Aj — \j)L/AE is 
used. Again, since iV'eoP = 0{e^), we have to order e 


P{ve v^) 


S 23 sin^ 2013 + 4e Jr cos 5 


(A+ - A-) - (Amgen -g) 
(A+ - Ao) 


Am? 


A 4 - - A_ 


sm 


2 ~ A-)A 


“h Sc Jr 


(Am; 


2 a3 

ren/ 


(A+-A_)(A+-Ao)(A_-Ao) 


. (A+-A_)L . (A_-Ao)A 
sm---sm---cos I 0 — 


AE 


AE 


AE 

I*- 

(3.14) 


here Jr, the reduced Jarlskog factor, is 


Jr = Ci2Si2C23S23Ci3S13- 


(3.15) 


This expression for the Ve —^ appearance channel probability is quite compact, despite 
that it contains all-order contributions of S 13 and a. In particular, it keeps the similar 
structure as the one derived by the Cervera et al. [7], which retains terms of order but 
is expanded by S13 only up to second order. 

This method of computing P{Uf. —)• v^) in the above offers the shortest path to the 
expression of the oscillation probability which is manifestly free from the apparent singu¬ 
larity as A_ —)• Ao because 1/(A_ — Aq) always appears adjacent to sin [(A_ — \o)L/AE]. 
We will refer this method as the “shortcut method” in the rest of this paper. 

The expression P{ve t'/x), eq. (3.14), is not quite of the form of eq. (2.1) because of 
the cos term in the expansion of cos ^<5 — ■ Tliis can be easily remedied, 

by using the following identity 

2sinA^_ sinA_ocosA+o = sin^ A+o — sin^ A_o — sin^ A^_ (3.16) 

which leads to exactly the form of eq. (2.1). However, if this is used then some of the 
terms are singular when A_ = Aq, yet the total expression is equivalent to eq. (3.14) and 
is finite. This is the reason why we prefer the form of eq. (3.14). 

The oscillation probability for i/g W channel can be obtained by all of the following 
three methods: ( 1 ) the similar calculation by the shortcut method, ( 2 ) using the unitarity 
relation P{i'e w) = 1 ~ P{^e ^e) — J’(^'e ^At)) (3) using the relation P{i'e w) = 

P{i'e —^ '■ C23 —>■ —'S23,'S23 —C23) [10], whose derivation is sketched in section 4.1. We 
note that Jr changes sign by the transformation. In general the oscillation probability in 
arbitrary channel can be obtained by ( 1 ) the shortcut method, or by ( 2 ) rewriting the 
expressions in appendix B using the formulas in appendix C. 


3.3.3 The general form of the Va —)• oscillation probabilities 

The expressions of the oscillation probabilities in the sector can be derived by one 

of the methods mentioned in the previous subsection. Then, one observes a remarkable 
feature that all the oscillation probabilities P{i'a (including the sector) can be 
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Ve Ve 

Ve V^ 

Vfj, Ve 

Ve Vr 

Vt Ve 

V^ Vr 

Vt V^ 

Vfj, V^ 

Ur 

Order e*^: 







Af_ 

-1 

sin^ 023 

COS^ 023 

- sin^ 023 COS^ 023 

- sin^ 023 

- cos'^ 023 

.aP _ .aP 
^+0 ~ ^-0 

0 

0 

0 

sin^ 023 cos^ 023 

- sin^ 023 COS^ 023 

- sin^ 023 COS^ 023 

Order ecosd: 







Bit = 

0 

1 

-1 

- cos 2023 

cos 2023 ~ 1 

cos 2023 + 1 

„aP _ „ap 

-“+0 ~ -°-o 

0 

0 

0 

- cos 2023 

cos 2023 

cos 2023 

Order esinA; 







gap 

0 

±1 

Tl 

±1 

0 

0 


Table 1. The values for the 5 coefficients for all oscillation channels, i^a —>■ and i^a —>■ J ^/3 to be 

used in conjunction with eq. (3.17). Note that they are 0, ± 1 or simple functions of 023- 


written in a universal form: 


P{Va Vp) = Sap 


{Af-} 


Z Z I 

Sa,Ca, + e 


fj^ap^ f J {(A+ - A_) - (Am2g„ - a)} 

{sy } (J„ cos J)- (A^ - a_)^(A^ - A„) - 


Sin 


2 (A + — X-)L 
4:E 


+ 4 


+ 4 


{^+o} + e {-S+o} (olrCOSii/cig) 

{^“ol 4 +e {-B-o} (Jrcosii/c^a) 


^ntren {(A+ " A-) - (Amge„ + o) } 
(A+ - A_)(A+ - Ao) 

Amgen {(A+ - A-) + (Amge„ + a)} 
(A+ — A_)(A_ — Ao) 


sin^ 


(A+ — Xo)L 
AE 


sin^ 


(A- - Ao)L 
AE 


-\- 86 Jrp 


(Amgen)^ 


(A+-A_)(A+-Ao)(A_-Ao) 


Sin • 


(A+ -A_)L 
AE 



— Xo)L 
~AE ' 


{C“^} cos^ cos^^^—{5“^} sinJ sin 

J AE ^ ^ AE 


, (3.17) 


where the eight coefficients , (7“^ and are given in Table 1. Notice that they 

are 0, ± 1, or the simple functions of 023 - 

In Table 1, we observe that three relations hold between the coefficients 


= A^q, 5“^ = Bfo and Bf_ = (3.18) 


These identities hold due to the invariance of the oscillation probabilities under the follow¬ 
ing transformation 

(g) —)• 7r/2-h (g) and A+-f-)-A_. (3.19) 

This invariance must be respected by the oscillation probabilities because the two cases in 
eq. (3.19) are both equally valid ways of diagonalizing the zeroth-order Hamiltonian. Look 
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at eq. (3.5) to confirm that it is invariant under (3.19). Then, the former two identities in 
(3.18) trivially follow, but the last one requires use of the kinematic relationship ^ 

sin A_|_sin A+o cos A_o = sin A_|_sin A_o cos A+o + sin^ A^_ 

We note that 023 dependence in the oscillation probabilities is confined into the five 
independent coefficients and all the other parameters, Am^’s, 0i2, 0i3, <5, E and a, are 
contained in the remaining part of the probabilities, which takes the universal form for all 
the oscillation channels. 

The antineutrino oscillation probabilities P{9a can be easily obtained from the 

neutrino oscillation probabilities. One can show that by taking complex conjugate of the 
evolution equation for anti-neutrinos of energy E that it is equivalent to the evolution 
equation for neutrino with energy equal to —E. Therefore, T’(i^a ■ E) = P{va ^p ■ 

-E). 

Here, we give some critical observations regarding the formulas (3.17) and the associ¬ 
ated Table 1: 

• The oscillation probability P{i'e —^ t'e), to first order in e, is of a simple two-flavor 
form with the matter-modified mixing angle 0 i 3 and the mass squared difference 
given by Am^gj^. As far as we know, having single two-flavor form as the expression 
of P{i'e t'e) in matter at this order is the unique case among all the perturbative 
frameworks so far studied. 

• Formulating the perturbation theory of neutrino oscillation with automatic imple¬ 
mentation of the oscillation probability formulas (2.1), which is to be carried out 
in section 4, is the principal result of this paper. However, the formula (3.17) goes 
one step further beyond eq. (2.1) by displaying the channel-independent universal 
function of L/E, up to the coefficients A±, B±, etc. which depend only on 623 . 

• With the universal form of the probabilities, unitarity can be trivially checked just by 

adding up the appropriate columns of Table 1 with proper care of using the correct 
signs in the last row for the intrinsic CP-violating terms.® One can easily confirm it 
by adding, for example, the t'e —> —>■ ^'/x, and Vg —>■ columns to give zero in 

each row. 

• If one uses the identity (3.16) to transform the sin sin cos 

term to the right-hand side of (3.16), then the transformational invariance given by 
eq. (3.19) will be manifest and eq. (3.17) will become the form of eq. (2.1). However, 
in the limit A_ —)• Aq some terms will diverge even thought the total expression is 
perfectly finite. 

^Notice that the relation 13“^ = C“^ needs to be satisfied only to order because these terms are 
already suppressed by e. 

® Here, we assume that the oscillation probability formulas are derived without using unitarity, excluding 
one of the options we mentioned earlier to derive them. 
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• In the universal form, the amplitude of a given oscillation frequency, to leading 
(zeroth) order in e, is controlled by </>. The amplitudes of the frequencies (A+ — 
Ao), (A_ — Ao) and (A+ — A_) are given proportional to Aj^qcos^ (p, A-osin^c/* and 

_sin^ ^ cos^ (/), respectively. For the channels involving Vf,-, only the (A+ — A_) 

frequency appears whereas in the — Vr sector all frequencies appear and are of 
equal magnitude near the ^13 resonance, when cj) ~ vr/4, but every where else either 
the frequency (A+ — Ao) ~ or (A_ — Ao) ~ Am^^^ dominates. 

Finally, we give more clarifying remarks on generic features of the oscillation proba¬ 
bility formulas presented in this section 3. 


• The term proportional to sin<5 in P{i'a —> z^/?), in eq. (3.17), is equal to 


8 < Jr sin 5 


(Am: 


2 \3 

ren/ 


(A+-A_)(A+-Ao)(A_-Ao) 


. (A+-A_)J . (A_-Ao)T . (A+-Ao)T 

sm-—-sm-—-sm 


AE 


AE 


AE 


up to an overall sign. This is because the quantity in {• • • } is just the Jarlskog factor 
in matter due to the Naumov-Harrison-Scott identity, [15, 16]. 


• Due to the above theorem sin <5 terms are always accompanied by the Jarlskog factor 
Jr = ci 2 Si 2 C 23 S 23 C^ 3 Si 3 . Though less well known, it can be proven that cos<5 terms 
must be proportional to the reduced Jarlskog factor Jrjc^z [17]. It can also be shown 
that both of the J-dependent terms must come with the suppression factor of e [15- 
17]. It is nice to observe that these general properties are realized in our formulas. 

• Association of the factor e and the all angle factors to the J-dependent terms indicates 
that they are genuine three flavor effects. In our general formula (3.17), all the terms 
with explicit e factor are of this J-dependent terms. Notice that other corrections of 
order e, though they do not show up explicitly, are implicitly contained in Am^en) 
A±^o and cj) which do not contain any 5 or O 23 dependence. 

• In vacuum, a = 0, the above oscillation probabilities reproduce the standard results, 
to first order in e. The form is somewhat unusual but we have checked that the 
expressions are identical. Also, we have checked that in the limit L/E ^ 0, where 
matter effects are negligible, the coefficients of the (L/E)"', for n < 3, terms in eq. 
(3.17), are identical to those in vacuum for all oscillation channels. 


3.3.4 Range of applicability of the formulas 

To discuss the range of applicability of our expressions, it is useful to first consider the 
vacuum expressions to first order in the expansion parameter e. For all channels, the 
expansion of the vacuum oscillation probabilities to first order in e does not include terms 
proportional to sin^ A 21 which starts at second order in e, 

sin2 A 21 = e^A^^^ + O(e^) ~ e^A^^ + O(e^). (3.20) 
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PC^e -> ^e) = 95,90,80.20,10,5% 





Figure 1. The iso-probability contours for the exact (solid blue) and approximate (dashed red) 
oscillation probabilities for upper left, Ve —> I'e, upper right, Ue —> and lower, —)■ v^. The 
upper (lower) half plane is for normal ordering (inverted ordering), whereas positive (negative) L/E 
is for neutrinos (antineutrinos). For treatment of antineutrinos, see section 3.3.3. The order of 
the contours given in the title is determined from the line L/E=0. The discontinued as one crosses 
Yep\E\ = 0 is because we are switching mass orderings at this point. In most of parameter space the 
approximate and exact contours sit on top of one another so the lines appear to alternate blue-red 
dashed. Note that, for L/E >1000 km/GeV and \ Yf , pE \ < 5 g.cm“^.GeV, the difference between 
the exact and approximate contours becomes noticeable at least for i/f. —>■ i/g and —>■ 

where Aji = Amj^L/AE. When, A 31 = 7 r/ 2 , that is at the first atmospheric oscillation 
maximum, ~ 0.002 which is small for the channels Vf, —?• Ux where x = e, or r 

since 1 — P{ve —t Ve), P{^e —t t'^) and P{i'e —t Vr) are all of order sin^ 2^13 « 0.1. However, 
at the second atmospheric oscillation maximum, Ag;^ = 37 r /2 and e^Agg ~ 0.02, which is 
significant compared to the sin^ 2^13 term. So in vacuum our first order expansion is only 
a good approximation for A|g < vr or L/E < 1000 km/GeV for these Ve channels. For 
the other channels, ^ u-r ^ Vt and Ufj, —>■ Ut, our hrst order expansion is a good 
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approximation to somewhat beyond L/E = 1000 km/GeV because the leading terms are 
not suppressed by the smallness of sin^ 26*13 • 

Then, what about the validity in matter? In section 5 we will argue that our per¬ 
turbative description is valid outside the solar resonance. Notice that the region without 
validity (no guarantee for approximation being good) is rather wide and includes the vac¬ 
uum because the solar resonance width |Aa| = \/ 3 (sin 20 i 2 /cos^ 0 i 3 )Am|^ is larger than 
the solar resonance position a = (cos 2 ^ 12 /cos^ 0 i 3 )Am|]^. We expect then that our helio¬ 
perturbation theory works for the matter potential a larger than a few tenth of | \. 

To give the reader a sense of the precision of our approximation we have plotted in 
Fig. 1, the contours of equal probability for the exact and the approximate solutions for 
the channels r'e —)• I'e —)• r'e and —)• The right (left) half plane of each panel of 

Fig. 1 corresponds to the neutrino (anti-neutrino) channel. As expected, for large values 
of the matter potential, \a\ > ^lAm^g^l we hnd we have no restrictions on L/E, to have 
a good approximation to the exact numerical solutions. Whereas for small values of the 
matter potential, |a| < ||Am^gj^| we still need the restriction LjE < 1000 km/GeV.® 

We note that most of the settings for the ongoing and the proposed experiments, 
except possibly for the one which utilizes the second oscillation maximum, fall into the 
region L/E < 1000 km/GeV. To improve the accuracy to larger values of L/E, especially 
for values of |a| < glAm^enl) second order perturbation theory in e is needed, which will 
be the subject of a future publication. 

3.3.5 An example of the power of our oscillation probabilities 

In this section we present an example of the power of our compact expression of the 
oscillation probabilities in understanding the physics of oscillations to first order in our 
expansion parameter e. For simplicity we consider the i'e disappearance channel. In fig. 2 
the I'e disappearance probability is shown, as a function of neutrino energy E, for baselines, 
L, of 3000 km (upper panel) and 5000 km (lower panel). As in hgure 1, the solid blue curves 
are drawn by using the exact expression of P{i'e i'e), and the dashed red curves by our 
compact formula in eq. (3.11). The black dotted curves are for the vacuum case. 

We first notice that our approximate formula agrees well with the exact expression in 
particular at higher energy. Secondly, because of the long baselines, the matter effect is 
sizeable, producing not only a large shift in the position of the hrst oscillation minimum but 
also an significant increase in the depth of the minimum. These two effects are correlated 
by the energy dependent quantity, (A+ — A_) which can be seen graphically in Fig. 3: the 
depth of the oscillation minimum changes from 

sin 2 2 % ^ ®”'20i3 

whereas the energy at which the the first oscillation minimum occurs changes from 

Am^enL (A+ — A_)L 

27r 27r 

® Of course, the boundary between these two regions should be interpreted as an approximate one. In 
fact, an exact boundary would sensitively depend on the definition of the difference allowed between the 
exact and the approximate probabilities. 
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Figure 2. The z/g disappearance probability as a function of energy E for baselines of 3000 km 
(upper panel) and 5000 km (lower panel). We have used the earth matter density 2.8 g/cm^. 

Because of the simplicity of our expression for P{ve eq. (3.11), these shifts are 

accurate to first order in the expansion parameter e. This simple understanding of the 
features of P(t'e —t i^e) is new to this paper. 

3.3.6 Comparison with the existing perturbative frameworks 

As we emphasized in section 2, our machinery has an advantage over the existing perturba¬ 
tive frameworks by having the minimum number of terms composed of sin [(Aj — \i)L/AE\ 
= 1,2,3) in the oscillation probabilities. This contrasts with the features of the exist¬ 
ing perturbative frameworks in which much larger number of terms than those minimally 
necessary as in (2.1) are produced. They include, typically, the terms with either extra 
LjE dependences or different frequencies in the sine functions, or often both, which easily 
obscures the physical interpretation. 

To give a feeling to the readers on how simple and compact our formulas are, we com¬ 
pare our expressions to the ones in the existing literatures to the same order in expansion. 
For definiteness, we pick the ones in ref. [10] to make the comparison, the most recent one 
among the reference list given in section 1. 

Our expression of P{ve —)• Ve) in (3-11) which has only a single term (ignoring unity) 
may be compared with eqs. (4.6) and (4.7) which consist of total 3 terms. With regard to 
P{ve if we count numbers of terms with different L/E dependence we have one 5- 

independent and 2 (5-dependent terms in (3.14), total 3. Whereas eqs. (4.8) and (4.9) in [10] 
have total 7 terms, 3 (5-independent and 4 (5-dependent ones. Our expression of Piy^j, —)• Vt) 
in (3.17) has 3 (5-independent and 2 (5-dependent terms, adding up to total 5. On the other 
hand, eqs. (4.10) and (4.11) in [10] contain total 8 5-independent and 10 5-dependent terms 
(2 sin5 and 8 cos5 terms), which adds up to 18. So not only do our expressions for the 
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oscillation probabilities have less than half the number of LjE structures, but the form of 
the L/E dependence is made manifest and identical to the vacuum form. 

We note that when our formulas are expanded by they agree, of course, 

with the existing ones calculated previously. Therefore, our formalism may be regarded 
as a systematic way of organizing the equivalent expressions to order into 

neater, structure-revealing ones. The benefit for having such simple expressions is that the 
physical interpretation of the formulas is transparent, as emphasized in section 2. 


4 Formulating the renormalized helio-perturbation theory 


In this section, we formulate the helio-to-terrestrial ratio perturbation theory, for short the 
helio-perturbation theory, which has the unique expansion parameter 


e 


Am2en ■ 


(4.1) 


We will show that use of its renormalized version is the key to the very simple formulas 
of the oscillation probabilities exhibited in section 3.3 and appendix B. In fact, there are 
two ways of deriving the oscillation probabilities, the S matrix method and the wave 
function method. Here we sketch both of them, leaving technical or computational parts 
into Appendices A and B. The meaning of the agreement between results obtained by both 
the S matrix and the wave function methods will be discussed at the end of this section. 

The S matrix describes neutrino flavor changes —)• Ua after traversing a distance L, 


t'o(L) = 5q,/3J//3(0), 

and the oscillation probability is given by 

l^a'i L) = 


(4.2) 


(4.3) 


When the neutrino evolution is governed by the Schrodinger equation, = Hv, S 

matrix is given as 


S 


Texp 


—i / dxH{x) 

Jo 


(4.4) 


where T symbol indicates the “time ordering” (in fact “space ordering” here). In the 
standard three-flavor neutrinos, Hamiltonian is given by 


/ 

1 

o 

o 

o 

_1 


a(x) 0 0 

' 

u 

0 Am^i 0 

+ 

0 0 0 


K 

0 0 Am|j^ 


0 0 0 

> 


(4.5) 


where the symbols are defined in an earlier footnote. For the case of constant matter 
density, the right-hand side of (4.4) may be written as We recapitulate here the 

earlier footnote: In (4.5) Am^- = — m? where m, denotes the mass of f-th mass 

eigenstate neutrinos. Position dependent function a{x) = 2\/2GFNf,{x)E is a coefficient 
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for measuring the matter effect on neutrinos propagating in medium of electron number 
density Ne{x) [1] where Gp is the Fermi constant and E is the neutrino energy. 

The neutrino flavor mixing matrix U is usually taken to be the standard form Upua 
given by Particle Data Group. We, however, prefer to work in a slightly different basis, 
for this paper, in which the flavor mixing matrix has a form (with the obvious notations 
Sij = sinOij etc. and 5 being the CP violating phase) 



’1 0 

0 


1- 

0 

0 

T —1 

1_ 

u = 

0 1 

0 

GpDG 

0 10 


0 0 


0 0 


10 0 ’ 


Cl 3 0 S13 

0 C23 S 23 p ^ 


0 10 

0 -5236“*'^ C23 


-Sl3 0 Ci3 


Cl2 Si 2 0 
-•512 C12 0 
0 0 1 


= U 23 U 13 U 12 


(4.6) 


under the understanding that the left phase matrix in the first line in eq. (4.6) is to 
be absorbed into the Vr neutrino wave functions. Being connected by the phase matrix 
rotation, it is obvious that our U and C/pdg give rise to the same neutrino oscillation 
probabilities. 


4.1 Choosing the basis for the renormalized helio-pertnrbation theory 

It is convenient to work with the tilde basis defined as Ua = (^ 23 ) 0 / 31 ^/ 3 , in which the 
Hamiltonian is related to the flavor basis one as 


H = uI^HU23, 


(4.7) 


where U 23 is defined in eq. (4.6). The S matrix in the flavor basis is related to the S matrix 
in the tilde basis S as 


S{L) = U23S{L)u) 


23’ 


S{L) = Texp 


—i / dxH{x) 

Jo 


(4.8) 


Notice that the matter term in the Hamiltonian (4.5) is invariant under the U 23 rotation. 

Hence, the dynamics of neutrino propagation in matter is governed by only the two mixing 

angles 612 and ^13 [20], which is also independent of 6 thanks to our convention of U in 

(4.6). From the first equation in (4.8), the relationships P{i'e —^ w) = P{i^e '■ C 23 —>■ 

—S 23 , S 23 —> C 23 ) etc. simply follow as used in section 3.3. 

The simplest formulation of helio-perturbative treatment in the tilde-basis includes 

decomposition of H into the zeroth and the first order terms in the expansion parameter 

Am2, 

Aid: as 


H{x) = 


2E 


a(x) I 2 n 
^ + ds 0 C13S13 


0 

Cl3'Sl3 


0 

0 


0 


U3 J 


+ 


Am? 


21 


Amii 


..2 „2 


'512^13 Cl2'Sl2Cl3 —'S12C13S13 


C 12 S 12 C 13 


G 2 


-C 12 S 12 S 13 


-S12C13S13 —Ci 2S12'S13 •S32S13 


(4.9) 
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To derive the compact formulas of oscillation probabilities we use slightly different 
basis, a renormalized basis, to formulate the perturbation theory. That is, we absorb a 
certain order e terms into the zeroth order Hamiltonian, H{x) = Hq{x) + Hi(x): 


Hq{x) 


Amren 

2E 




“(a;) 1 „2 

0 

Cl3'Sl3 


r 

*12 

0 

0 ’ 


0 

0 

0 

+ e 

0 

^12 

0 


C 13 S 13 

0 

^13 J 


0 

0 

*12 _ 

> 


Hi{x) = eci 2 Si 2 


Am^-en 

2E 


0 Ci3 0 

Cl3 0 -S13 

0 -Si3 0 


(4.10) 


(4.11) 


where Am^g^ = Am|]^ — ^ defined in (3.2) and (3.3). H{x) 

with (4.10) and (4.11) is identical with the tilde-Hamiltonian in (4.9). Note the simplicity 
of the perturbing Hamiltonian, Hi, especially the zeros on the diagonal elements. The 
diagonalization of Hq leads to the eigenvalues given in section 3.1, eq. (3.1), and because 
of the zeros along the diagonal of Hi, there are no corrections to these eigenvalues at hrst 
order in a perturbative expansion. 

The authors of [5] treat the order e effect in the Hamiltonian as a renormalization of 
the matter potential, whereas we regard it as a renormalization of Aml;^. 

Our treatment can be easily generalized to cases with matter density variation as far 
as the adiabatic approximation holds. See e.g., [10] for such a treatment. However, we 
prefer to remain, for simplicity and compactness of the expressions, into the formulas with 
constant matter density approximation in the rest of this paper. 


4.2 Hat basis 

We transform the Hamiltonian H = Hq + Hi, from the “tilde” basis to the “hat” basis, 
using 


Ho = uIHoU^, Hi = 


where the unperturbed Hamiltonian Hq is diagonal, 


Ho 


1 


A_ 0 0 
0 Ao 0 
0 0 A+ 


(4.12) 


(4.13) 


We take the following form of unitary matrix to diagonalize Ho'. 

cos (j) 0 sin (j) 


U^ = 


0 10 


(4.14) 


— sin (p 0 cos (p 

The expressions of the zeroth order eigenvalues A_, Aq, and A+, are given in eqn (3.1). 
Similarly, cosine and sine p are given in eqn. (3.5). 
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Also the perturbing Hamiltonian, Hi^ retains it’s simple form thanks to that the 
rotation keeps “zeros” in Hi unchanged both on the diagonal and in the top-right and 
bottom-left corners, 


Hi = uIHiU^ 


= eCi2Sl2 


2E 


0 cos ((?!) — 6 * 13 ) 0 

cos ((/> - ^is) 0 sin {(j) - 013 ) 
0 sin {(p — 013 ) 0 


(4.15) 


In fact, Hi is identical to Hi with 0 i 3 replaced by ( 0 i 3 — (p). 

4.3 S Matrix and the oscillation probability 

The S matrix in the flavor basis is related to the S matrix in the tilde and the hat bases as 


S{L) = U23S{L)uI^ = U23U^S{L)UlUl, 


Ttut 


(4.16) 


where we have used explicitly the fact that the matter density is constant: 


S{L) = Texp 


—i / dxH{x) 


= U^T exp 


—i / dxH{x) 

Jo 


UI^U^S{L)UI. (4.17) 


To calculate S{L) we define Q{L) as 

n{L) = 

Then, n(L) obeys the evolution equation 

i—Q(x) = Hiid(x) 
dx 

where 


(4.18) 


(4.19) 


Hi = 


(4.20) 


Then, 17(x) can be computed perturbatively as 

n{L) = l + {-i) [ dxHi{x) + 0{e‘^). (4.21) 

Jo 

Collecting the formulas the S matrix can be written as 

S{L) = t/23t/0e-*^°^17(L)t/Jt/j3 (4.22) 

Thus, we are left with perturbative computation of 17(L) with use of (4.20) to calculate 
the S matrix. With the S matrix in hand it is straightforward to compute the oscillation 
probabilities by using (4.3). We leave these tasks to appendices A and B. 
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4.4 Mass eigenstate in matter: V matrix method 

In this section we calculate the I/-matrix directly using our perturbation theory. If we 
switch off the perturbation Hi, the mass eigenstates in matter, to lowest order, are given 
by the hat-basis wave function which are the eigenstates of Hq in (4.12), and since 
Hq is diagonal, we have 

^ (4.23) 

Thus, the V matrix is given to zeroth order by = ?723t^0 whose explicit form is given 
in section 5.3, and also in eq. (3.8). 

In order to obtain the mass eigenstates in matter to first order in e, Ui = let 

us compute the first order correction to the hat basis wave functions. Using the familiar 
perturbative formula for the perturbed wave functions 

,(l) ^ ^ pf (4.24) 

... At 

with Hi in (4.12), and the Aj’s are given by the eigenvalues of Hq, see (3.1). Then the 
mass eigenstate in matter z/j can be written to first order in e as: 


("-] 

( 

\^+) 

d 


1 

Cl 2 
0 


^^''hen A_-Ao 


,A™2 ci2Si2C(^_e^3) 

A_-Ao 


e^ZZtren A+-Ao 


0 

C12 

1 




( 0 ) 


e^ZZtren A+-Ao 


V 


>r 

40) 


Using (4.23), this equation is of the form Vi = V^Ua which can be inverted to easily obtained 
the U-matrix given in eq. (3.7), 




' 

O 

T—1 

1 

o 

o 

o 

o 

1 

U = U23U^ 

1 + eci2Si2Am^gj^ < 

A--Ao 

10 0+ 

1 — An 

0 0 1 




< 

c 

c 

c 

c 

1 

c: 

J 


(4.25) 


This can be used to directly compute the oscillation probabilities as was performed in 
section (3.3), or by inserting the V matrix elements into the expressions of the oscillation 
probabilities in (4.6). This V matrix method was used to calculate the matter effect 
correction in the oscillation probabilities [18]. For a different approach toward simpler 
expressions of the V matrix elements see [19]. 

In sections 4.3 and 4.4, we have sketched the two different methods for calculating 
the oscillation probabilities, the S'-matrix method (section 4.3) and the U-matrix method 
(section 4.4). The results obtained by the both methods agree with each other, and the 
expressions of the oscillation probabilities are of the form given in eq. (2.1). 

Let us make a comment on what happens if we use a different decomposition of the 
Hamiltonian into unperturbed and perturbed parts as in (4.9). Then, the perturbed Hamil¬ 
tonian has nonzero diagonal terms, and this produces first order corrections to the eigen¬ 
values. If we use these expressions for the eigenvalues and expand by the small expan¬ 
sion parameter we obtain terms that do not exist in (2.1), such as the ones proportional 
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Figure 3. The flow of the three eigenvalues in the normal mass ordering, abbreviated as NO, (left 
panel) and in the inverted mass ordering, abbreviated as 10, (right panel). The exact eigenvalues 
are depicted with colored lines, green for vi, red for V 2 , and flight blue for v^. The eigenvalues 
calculated by using our renormalized helio-perturbation theory are drawn by the black dashed lines 
whose state labels are marked on the figures. The approximate eigenvalues for z/q and v- cross at 
the solar resonance, whereas the exact eigenvalues for vi and V 2 repelled at solar resonance. To 
make these features visible we used a solar Am|]^ three times as large as the measured value. 

to sin [(Aj — \i)L/2E]. Of course, the S-matrix method will give the same expressions. 
Therefore, absence of the explicit first-order correction to the eigenvalues is essential to 
keep perturbative expressions of the oscillation probabilities of the form (2.1) to order e. 
In our case it is guaranteed by vanishing diagonal elements of the perturbed Hamiltonian 
(4.11). 

5 More about the renormalized helio-perturbation theory 

In this section, we critically examine the framework of the renormalized helio-perturbation 
theory. Despite a drawback of the current framework (which is to be described below) we 
argue that our perturbation theory works apart from the vicinity of the solar resonance 
crossing. 

5.1 Exact versus zeroth-order eigenvalues in matter 

The three eigenvalues of the Hamiltonian are written as where i runs over 1,2,3 for 
the exact eigenvalues, and i = -,0,-1- for the zeroth-order eigenvalues in our perturbative 
framework. In Figure 3 the A* are plotted as a function of a, the Wolfenstein matter 
potential, for both the exact and the zeroth-order ones given by eq. (3.1). It is clear 
in Fig. 3 that our zeroth-order eigenvalues fail to treat the solar-Am^ scale level crossing 
correctly. As one goes through the solar resonance in the exact solution the two eigenvalues 
involved, the red and green, repel one another, whereas in our perturbative solution the 
two corresponding eigenvalues cross with each other. 

We point out here that the drawback is the common feature among the similar pertur¬ 
bative treatments of neutrino oscillation available in the market. In this section we argue. 
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for the first time, that despite the drawback, our perturbative framework successfully treat 
the flavor content of these two states in region reasonably far from the solar resonance. 

We first note that, despite the feature, the atmospheric and solar resonances occur at 
the correct values of the matter potential with our zeroth-order eigenvalues. The atmo¬ 
spheric resonance occurs when A+ — A_ is a minimum, which is at 

a = cos 2613 

and the minimum difference between A+ and A_ is Am^gj^ sin 2^13 as expected. The solar 
resonance occurs when A_ — Aq is a minimum; this occurs when 

a = eAm^gn cos 2612 / cos^ ^13 

and the minimum difference is zero! This is not the value determined by the full Hamil¬ 
tonian which is ~ eAm^en sin 2012. Therefore, while our perturbative scheme treats the 
atmospheric resonance correctly to order e, it misses the effects of the solar resonance. 

5.2 Eigenvalues in vacuum and in the asymptotic regions a —)• ±00 

In this and the next subsections, we will give the arguments to indicate that our renormal¬ 
ized helio-perturbation theory works apart from the vicinity of the solar resonance despite 
the issue mentioned above. 

We first show that the zeroth-order eigenvalues given in ( 3 . 1 ) agrees with the exact 
ones to order e in the asymptotic regions a —)• ±00. We use the characteristic equation of 
the full Hamiltonian ( 4 . 5 ) to derive (i = 1 , 2 , 3 ) 

Aj = (a -|- Am|]^ -|- Am^i) , 
i 

^^XiXj = Am^iAmli-|-a {(ci2-|-Si2Si3)Am2i-|-Ci3Am|i} , 

A1A2A3 = c^2Ci3«Am2iAm|3. ( 5 . 1 ) 

Then, by using the asymptotic expansion of A’s one can obtain to leading order in the 1 /a 
expansion (at a —)• -|-cx) in the NO case) 

A2 = cf3Am§i -h sf2Sl3^"i21 = (c?3 + 

Ai = cf2Amii = c?2eAmrg„, 

A3 = a -h Si3Am|i -k sf2Cl3^"^-21 = “ + (■sL + Sl2e)^"l?en- ( 5 - 2 ) 

At a —)• —00, A+ and A_ must be interchanged. They are identical to the ones computed 
with our zeroth-order eigenvalues given in ( 3 . 1 ): 


(x.\ 

NO 

/ 

(c?3 + Si2e)Am2gn \ 

Xo 

= 


c?2eAm2gjj 


a=+oo 


+ ('S?3 + •Sl2e)'^"lren / 
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(5.4) 



NO 

'^« + ('S?3 + 'Si2e)Am2g^\ 

Ao 

= 



a =—00 

\ (c?3 + Si2e)Am2gn y 


For the 10, the analytic expressions of Aj {i = —, 0, +) at a —)• ±cx) is the same as those of 
Aj at a —)• =Foo for the case of NO.^° 

When the matter potential vanishes, a = 0, the mass squared eigenvalues are given for 
the both mass orderings by 



( s\^eAml^^ \ 

Ao 

c?2cAm2g„ 

\^+/ a=0 

\ / 


whose order e terms are different from the exact values, Ai = 0, A 2 = and A 3 = 

To understand the meaning of the failure at order e we need to discuss what happens 
to the flavor content of the matter mass eigenstates as the matter potential changes from 
below to above the solar resonance. This will be done in the next subsection. 

Similarly, the asymptotic behaviour of the angle </>, i.e., @13 in matter can be easily 
worked out. With the definition of (p in (3.5), it is easy to show that 4> takes on the following 
values as a is varied from — 00 to + 00 . In the case of normal mass ordering, 


4>no 


0 , 

^13, 


7 r 

4’ 



ZE 
2 ’ 


a = — 00 
a = 0 

a = cos 26*13 

a = 2 Am 2 gn cos 2013 
o = + 00 . 


while in the case of inverted mass ordering. 


Pio 


TT n 

2 ^13? 

TT 

4 ’ 
^13 
0 , 


a = — 00 

a = 2 Amren cos 20 i 3 
a = Am^gjj cos 20 i 3 
a = 0 
a = + 00 . 


(5.6) 


(5.7) 


It reflects a natural view that physics at a —)• +00 for the normal mass ordering corresponds 
to the one at a —)• —00 for the inverted mass ordering at least to leading order in e. 


5.3 Neutrino mixing matrix in matter and the v flavor content at o —?• ±00 

Since the two levels cross at the solar resonance in our perturbative treatment, one may 
expect that our treatment fails completely beyond the solar resonance, i.e., in the region 
with matter density higher than the resonance. However, we will show in this subsection 
that the flavor contents of the three eigenstates are correctly reproduced at least in the 

In fact, the first line in (5.1) holds exactly. 
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asymptotic region. That is, the zeroth-order V matrix describes correctly the asymptotic 
behaviour of the exact eigenstates in matter. 

Suppose we denote the exact flavor mixing matrix in matter as the (almost standard) 
U matrix defined in (4.6). Let us discuss the case of NO first. At a —?• —oo U is nothing 
but in (3.8). Notice that S 23 in matter is frozen to its vacuum value for e < 0.1, and 
S 12 ~ 0 at a —)• —00 [3]. At o —)■ + 00 , S 12 — 1 and C 12 — 0 the matter U matrix is identical 
to (3.8) if we interchange vi and 1^2 apart from re-phasing factor —1 for the new second 
mass eigenstate. 

To make the meaning of this feature clearer we write down here the flavor content of 
the states Ui {i = 1,2,3) at o —)> ± 00 . Noticing that Ui = the flavor composition 

at a —>■ —00 is given at zeroth order by 

= C^Ve - S^{s23e~''^l'f, + C23^'r) = 1^- 
= {c23^^^l - S23e*^^'T} = I'O 

z^3 = S^Ve + C</,(s23e“*‘^t'/, + C23^'r) = l^+ (5.8) 

Whereas the composition at a —)• +oo is given by 

= -{c23^^l - •5236*'^ 

= C^Ve - S^{s23e~'‘^l'f, + C23W) = 1^- 

Z^3 = S^Ue + C^{s23e~'^^fi + C23W) = 1 ^+ (5.9) 

The flavor compositions given in (5.8) and (5.9) imply that the flavor content of the lower 
two mass eigenstate in matter is correctly described in our perturbative framework despite 
the failure of describing the solar level crossing. Note, the combinations of and Vr in 
the (• • •) and {• • • } are orthogonal in the above expressions. 

In the case of 10, essentially the same discussion goes through. The asymptotic be¬ 
havior of 012 at a —)• Too is the same as that in NO. The asymptotic behavior of 0 i 3 at 
o —>■ Too for NO is mapped into the one at a —>■ Too for 10. But, this is already taken care 
of by the definition of cf) given in (3.5). Therefore, the same U matrix in matter as the one 
in NO are obtained at a —)• Too. Then, the same flavor compositions as in (5.8) and (5.9) 
follow for 10. Again it is consistent with those we expect from the level crossing diagram. 

We add that our state always corresponds to 1/3 state. At a —)• Too in NO and at 
a —>■ —00 in 10 the electron neutrino component is all in this = 123 state. At a —)• —00 
in NO and at a —)• Too in 10 the electron neutrino component is all in {vi for NO, 1^2 
for 10) state. In vacuum for both NO and 10, = S 13 , which is the correct value 

for the V 3 = state. 

To summarize: Despite that the eigenvalues calculated by our helio-perturbation the¬ 
ory do not show the correct behavior at around the solar level crossing, the flavor composi¬ 
tion of the states are correctly represented by our zeroth-order states. It is the very reason 
why our perturbative formulas work at much higher and lower values of Y^pE compared 
to the one at the solar resonance. This point escaped detection in the previous literature 
to our knowledge. We have seen that they work practically whole regions except for the 
vicinity of the solar resonance. 
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6 Concluding remarks 


We have developed a new perturbative framework of neutrino oscillations, which allows us 
to derive compact expressions of the oscillation probabilities in matter, which we examined 
in this paper to order e = ~ Am^j^/Aml^. Although vacuum like in their 

forms, our compact results keep all-order effects of S 13 and the matter potential a, while 
using e as the unique expansion parameter. 

The characteristic features of our perturbative framework (dubbed as the renormalized 
helio-perturbation theory), in contrast to the ones previously studied by various authors, 
are as follows: 

• We use the renormalized basis (4.10) defined with the atmospheric mass-squared 

difference corrected by an order e quantity, Am^en = ~ decom¬ 

position of the Hamiltonian into the unperturbed (4.10) and perturbed terms (4.11), 
is done in such a way that there is no diagonal entries in the perturbed one. It makes 
the eigenvalues of the zeroth order Hamiltonian the correct ones to order e. Usage 
of this decomposition is crucial to obtain the form of the oscillation probabilities 
(2.1) akin to the form in vacuum. This facilitates an immediate physical interpreta¬ 
tion that the frequencies of the oscillations are determined by (A+ — A_), (A+ — Aq), 
and (A- — Aq), where the A’s are the eigenvalues of the zeroth order Hamiltonian, 
see eq. (3.1). They determine the oscillation pattern, the functional form oi L/E 
dependence, for all the oscillation probabilities as in eq. ( 2 . 1 ) to this order. 

• The amplitudes of the sin [(Aj — Xi)L/4:E] terms in the oscillation probabilities, (2.1), 
can be determined either by the iS-matrix method, or the U-matrix method. The 
latter may be simpler because of no first order correction to the eigenvalues. In 
contrast to the complexity of the exact solution, this perturbative approach ensures 
calculability of the V matrix elements in the form of power series, see eq. (4.25). It 
leads to the simple analytic expressions of the oscillation probabilities which have 
a universal form to first order in e, as shown in eq. (3.17). The coefficients in this 
universal form are 0 , ±1 or simple functions of 823 , see table 1 . 

• From the universal form of the oscillation probabilities, eq. (3.17), we readily observe 
the followings: 

(1) Each one of the zeroth-order terms in the oscillation probabilities, in all channels, 
takes the same form as the corresponding two-flavor oscillation probability in vacuum 
but with use of the eigenvalues and 613 in matter. In the special case of z^e 

only the term with the particular frequency oc A+ — A_ is non-zero, so the LjE 
dependence is of the two flavor form. The dominance of this frequency also occurs 
in the oscillation probability for z/g —)> and Ue ^ I'r- 

(2) All the first-order correction terms with the explicit factor e (~ AmQ/Am^) 
in the oscillation probabilities are proportional to either cos(i or sin (5. They are also 
suppressed by the angle factor oc C 12 S 12 C 23 S 23 S 13 . Hence, they are the genuine three- 
flavor effects. The explicit expressions of the 5-dependent terms in our formulas are 
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in agreement with the general theorems. 

(3) Unitarity can be trivially checked. 

Further comments are made in section 3.3.3. 

Despite the success of our current framework in almost all regions including the one 
around the first oscillation maximum, it has a clear drawback. It fails to accommodate the 
physics at around the solar-scale resonance. This is related to the unphysical feature that 
the two eigenvalues (A_ and Aq) cross with each other at the solar resonance, the universal 
fault in all the perturbative treatment which include the expansion parameter e in the 
market. However, we have offered, for the first time to our knowledge, an explanation why 
our framework works beyond the resonance despite the failure in treating the solar level 
crossing. We hope that we can return to the problem as a whole in the future. 

As we noticed at the end of section 3.1 our renormalized Am^en is identical, to order 
e, to the effective Am^ measurable in an (anti-) Vf. disappearance experiment in vacuum. 
It is a tantalizing question whether it is just a coincidence, or is an indication of something 
deep. 

It is quite possible that the features of the current framework mentioned above natu¬ 
rally generalizes to higher orders. That is, one can demand that the oscillation probabilities 
of the general form in ( 2 . 1 ) calculated by the H-matrix method, with the carefully chosen 
zeroth-order eigenvalues, agree with the ones from the ^-matrix method and are correct to 
certain order in e. Our result in this paper may be regarded as an existence proof of this 
concept to order e. Since we know that this is true in the exact form of the oscillation prob¬ 
abilities (assuming adiabaticity) [3], it is likely to be correct in each order in perturbation 
theory. It may or may not require higher order renormalization in Am^en- 


A Calculation of S matrix elements 
A.l Computation of S matrix elements 

By using il\ in (4.15) and = diag g*Aoa;/ 2 E^ giA+a;/ 2 £;j ^ one can easily com¬ 
pute Hi[x) = Then, using eq. (4.21), the first order term of D(x) can be 

calculated as 


Di(a:) = [—i) f dx'Hi{x') 

Jo 


= -ieAmrenCl2Sl2 


0 

g—i(A_-Xo)a;/2E_ 

H4>-Si3) -j(A_-Ao) ’ 

0 


gi(A_ -Ao)a:/2B_j^ 

i(A_-Ao) ’ 
0 


gi{A_|_-AQ)a:/2B_j^ 

^i4>-Sl3) i(A+-Ao) ’ 


0 

^i<f>-Sl3) -i(A+-Ao) 

0 


(A.l) 


where we have introduced the simplified notations, = cos ((/> — ^ 13 ), = 

sin ((/) — ^is), etc. The simplicity in the structure of (A.l) with many zeros is the math¬ 
ematical reason why the expressions of neutrino oscillation probabilities are so simple in 
our renormalized helio-perturbation theory. 
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The S matrix is given hy S = e = e [1 + f^i(a;)] where 1 denotes the unit 

matrix. Then, the elements of S matrix are given by: 

Q —iX—xf2E Q _ —iXqx j2E q _ —iX^x/2E 

^ee — ^ 1 ^ ’ *^rr — ^ 5 

Sex ~ Sq-e — 0, 

^^—iXQx/ 2 E _ ^—iX-x/ 2 E^ 


Sefi — Sfj^e — ieA?77-j.gj^Ci25l2C(0_0^3) ■ 

S^q- Sq-^ i6^77lj.Qy^Ci2^12^(^(p—0i3)~ 


i{X- - Ao) 

'^—iXoxl2E _ ^—iX^x/2E^ 


i(A+ - Ao) 

A.2 The relationships between S, S, and S matrices 

The relationships between S, S, and S matrices are summarized as 


S{x) = U^S{x)ul S{x) = U 2 sS{x)E 


23 - 


The left equation can be written explicitly as 


S = 


See Se^ Ser 
S^e S^^ S^r 
Sre SrfM Srr 



^(p 0 S(p 


= 

0 1 0 



1 

c« 

0 

1_ 



See Set, Ser 
SfMe S^^J, S^T 
Sre Sr^^ Srr 



^(p 0 ^(p 


0 1 0 


c« 

0 

Ci 

1 _ 


Hence, the elements of S can by written in terms of the elements of S as 
See — ^(jySee T ^(pSrr T ^Ser T Sj- 

Se^ — -|- S(pSq-^^ 

Ser — ^(jqSer ^(jySre C050 ^*§ee Sq- 
Sfj,e — ^(pSfj^e H“ ^(pS^q- — Sefj^ 

S^r = - SfpSfj^e T (^rpS^T, 

Sre — S^Ser T C^Sq-e CfpS/j) ^See S^r ) — ‘S'er; 

St^ — SfpSefj, T CfpS-j-fj, — SfiY 
Stt — S^See T C-^Stt (Ser T S^ 


(A.2) 


(A.3) 


(A.4) 


(A.5) 
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Similarly, the elements of S can be written in terms of the elements of S matrix as 



(A.6) 


Then, the oscillation probabilities are simply given by P{vi 3 —)• eq. (4.3). 

B Expressions of neutrino oscillation probabilities 

With the expressions of S matrix elements obtained in appendix A, and using eq. (4.3) 
it is straightforward to calculate the neutrino oscillation probabilities. Similarly, one can 


insert the V matrix elements given in (3.7) into (2.1) to obtain the equivalent results. In 


this appendix we only give the resulting expressions, but of all the oscillation channels 
under the viewpoint that unitarity is not to be imposed but must be proven to show the 
consistency of the calculation. To compare these with the results of section (3.3), use of 
the identities given in Appendix (C) may be of some help. 

The only comment worth to give here is about the simple method for transformation 
C 23 —^ —S 23 and S 23 —^ C 23 to obtain P{i'e w) from P(z^e —>• or Piy-j- —>• Ut) from 

P{Vfi —)• r'/i), which is utilized in section 3. Though we work with the rephased flavor 
mixing matrix dehned in (4.6), the transformations produces Ser from S'e^, and Srr from 
up to an overall phase, see eq. (A.6). 

B.l Oscillation probabilities in ly^—row 
P{i'e ^e) is extremely simple as 



(B.l) 
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where L is the baseline. The reasons for the simplicity is discussed in depth in section 3.3. 
P{ve —>• Vp) and Pi^e ^t) are given by 

P{ve i^p) = ^clslsls sin^ 


AE 


+ 4eAmrgnCi2Si2C23S23C</,S</,COs5 
1 






1 


sm 


2 i^+ ~ A-)T 


AE 


.2(A+-Ao)T .2(A_-Ao)T'1 

sm -—-sm - 


AE 


AE 
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— Ao) 


. (A+-A_)L . (A_-Ao)L . (Ao-A+)L 

X sm---sm-—-sm 


AE 
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(B.2) 


or . ^ ^ 2 2 2-2 ('^+ A_)L 

P{Ve l^r) = 4c^S^C23 S™ 


AE 


- 4eAmrenCl2Sl2C23S23C</,S</,COS(5 
1 


■C(/)C(^_033) - 


(A_-Ao)^"‘^"(‘^-'^^\a+-Ao) 

. 1 1 


1 


sm 


2 i^+ ~ A-)T 


AE 


■ 2 ('^+ ~ ^o)L . 2 ~ ^o)L \ 

sin - — -sm ---> 

AE AE I 


+ 8eAmrenCl2Sl2C23S23C</,S^sin J |c<^C( 0 _ 033 )-^^^— _ + S0S(^_e33)-^^^^—-^io) } 

. (A+— A_)L . (A_ — Ao)T . (Aq — A+)L m Q^ 

- AE -- AE -- AE -• 

It is almost trivial to verify unitarity in the Ue—Tow: P{ve ^e) + P{^e i^p) + P{i^e 

i/P = 1. 


B.2 Oscillation probabilities in Vp—row 

P{i'p —)• fe) is related to the T-conjugate channel probability P{i'e —t ^fi) as P{vp —>■ 
Ue 6) = P{i'e —t Vp : —(5), whose latter can be obtained by replacing 5 by —6 in (B.2). 
Therefore, we only give the expressions of P{i'p —)• Vp) and Pivp —>■ Vt)'- 

P{i^p Vp) 


= 1 - 4523^^0 


■ 2 (-^+ “ A_)L 


“ 4 C 23 S 23 


AE 
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(B.5) 


AE AE AE 

With the above results the unitarity in row can also be verified: —)• I'e) + 

v^) + P{v^ Ur) = 1. 

B.3 Oscillation probabilities in Ur—row 

P{vr —)■ Ve) and P{ut —t u^) can be given by their T-conjugate channels: P{ur Ua '■ 6 ) = 
P(ua -7- Ur : —(5). Therefore, we only give the expressions of P(ur —)• Ur) below. 


P{Ur Ur) 


,4 „2„2„- 2(A+ A_)L 2 „2 
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(B.6) 


The unitarity in row can also be verified: P{ur —)• Ug) + P{ur —)• u^) + P{ur —)• i^'t) = 1- 


C Some useful Identities 

A straightforward derivation of the general expressions of the oscillation probabilities con¬ 
tains the following form of CP- or T-violation terms 

- 2 ^ lm[Vc..V^,V:^V 0 j] sin , (C.l) 

j>i 
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where i = 0, —, + (or 1,2,3). To cast this term to the one in (2.1), one needs the following 
identity 


. (A+-A_)L . (A+-Ao)L , . (A_-Ao)L 

sin--sin--h sin ■ 


2E 


2E 


2E 


^ . (A+-A_)L . (A_-Ao)L . (Ao-A+)L 
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4£; 


AE 


(C.2) 


Here, we list some more formulas which may be useful to understand the relationship 
between different expressions of the oscillation probabilities: 
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